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We discuss analytical properties of partial waves derived from projection of a 4-legged amplitude 
with crossed-channel exchanges in the kinematic region of the direct channel that corresponds to the 
XYZ peaks in charmonium and bottomonium. We show that in general partial waves can develop 
anomalous branch points in the vicinity of the direct channel physical region. We show that this 
effect only occurs if masses of resonances in the crossed channel are in a specific, narrow range. 
We estimate the size of threshold enhancements originating from these anomalous singularities in 
reactions where the Z c (3900) and the Zb(10610) peaks have been observed. 


I. INTRODUCTION 

There is significant interest in the physics of heavy quarkonia stimulated by discoveries of narrow peaks in the 
spectrum. Such peaks may indicate existence of new hadrons. A recent review of the experimental situation and 
of the various theoretical models can be found, for example, in pQ. While the quark model provides a remarkably 
accurate description of the heavy quarkonium spectrum, these new peaks, also referred to as the XYZ states, appear 
at masses that do not, in a natural way, derive from the quark model [2]. Several of these peaks have been observed in 
invariant mass distributions of meson pairs that contain one heavy quarkonium, e.g. the J/tp or the T, and one light 
meson. The possibility that these new hadrons are therefore multi-quark bound states has been explored, for example 
in mm. The peaks appear near thresholds for production of meson pairs with open flavor, e.g. DD* or BB *, and for 
this reason it has also been suggested that binding between the two flavored mesons may be responsible for generating 
some of the XYZ’s [ME]. It should be recognized, however, that distinction between multi-quark and hadron bound 
states is complicated [9] and in any case requires sophisticated amplitude analysis and precise data m ■ It has also 
been suggested m and studied in HME] that coupling to nearby open channels may produce peaks even without 
presence of new hadrons. In the mathematical language of amplitude analysis the analogous statement would be that 
threshold cusp may be enhanced not only by a nearby pole but by a carefully arranged set of branch points, also 
referred to as anomalous thresholds. Amplitude poles correspond to physical particles (stable or unstable) while cuts 
represent the effect of ”forces”, i.e. exchanges of known particles in the crossed-channels. If the XYZ peaks were due 
to the latter it could potentially weaken the case for new hadrons. To address this issue it is therefore important to 
perform a systematic analysis of the ” forces” and the associated singularities. 


II. THE MODEL 

We consider a decay of a heavy object of mass M, hereafter referred to as the M particle, to a quasi stable heavy 
state of mass M' and a pair of degenerate, light mesons of mass p. We are interested in the energy dependence of the 
low spin partial waves, in particular the S'-waves, in the M'p channel. In this section we discuss a generic case and 
in the section that follows we consider threshold behavior in J/tttt and T(1S)tt production from Y(4260) —>• J/ip 7r7r 
jl9H2T] and Y(5 S) -A T(1S')7T7t [22 23 decays, respectively. 

The reaction of interest is M(p\ ) -A M'(ps) + p{p 2 ) + t{pa) with pi referring to the 4-vectors. Ignoring spin, the 
reaction amplitude, A(s, t ) is a scalar function of two independent Lorentz invariant Mandelstam variables, which we 
choose as, s = (p% — P 2) 2 and t = (p\ — P3) 2 . They correspond to the invariant mass squared of the M'(ps)p(p 4 ) pair 
and momentum transfer squared between the two heavy mesons, respectively. The amplitude in the kinematics of the 
decay region, s > ( M' + p) 2 , t > 4/i 2 can be obtained by analytical continuation in t of the amplitude describing the 
5-channel scattering process, M(pf) + /i(—P 2 ) M'(pz) + p{p 4 ). It is the t-channel singularities of the latter that 
are responsible for anomalous singularities of the 5-channel partial waves. Similar analysis applies to the ^-channel 
where u = (pi — P 4) 2 . Bose symmetry requires the amplitude to be 5 H w symmetric. 

The 5-channel partial waves, Ai(s), describe production of the M' p system in a state of fixed angular momentum, l. 
Ignoring spin of the external particles, partial waves are the coefficients in expansion of A(s,t) in a series of Legendre 
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polynomials, 


A(s,t) = I6n^2(2l + 1 )Ai(s)Pi(z 3 ). 
1=0 


( 1 ) 


The argument of the rotational functions, = z s (s,t), is the cosine of the scattering angle in the center of mass 
frame of the 8-channel. When considered as complex functions of the invariant mass squared, 8, the partial waves 
Ai(s) have branch points at all 8-channel production thresholds. The cuts associated with these branch points are 
located on the positive real s axis and unitarity determines the discontinuity across the cuts. These cuts are referred 
to as the right hand cuts. Through the partial wave projection, 


Ai(s) 



dz s Pi(z a )A(s,t), 


( 2 ) 


threshold singularities of A(8, t) in the t and the u channel also lead to branch points in the complex s plane. The 
associated cuts are generically referred to as the left hand cuts, even though, as it will be the case here, they don’t 
always he on the negative real 8-axis. In the following we consider only the S -waves (l = 0) and drop the angular 
momentum subscript on the partial waves. Higher partial waves are suppressed at thresholds due to the angular 
momentum barrier factors. 

In the cases considered here, y/s is of the order of several GeV and there are several open channels contributing 
to the right hand side discontinuity. As long as their thresholds are far from the region of interest they will not lead 
to rapid variations in 8. It is therefore sufficient to focus on the channels, which open near the XYZ production 
thresholds. For simplicity we assume that the relevant channel consists of two identical mesons with mass m. In the 
cases considered in the following section these will correspond to the BB* or the DD * state. Here, for simplicity, we 
ignore the mass difference between the pseudoscalar and vector mesons. We assume that masses satisfy the following 
relation, M > 2m and 2 m > M', i.e. we take mass of the decaying/produced quarkonium to be slightly above/below 
the open flavor threshold. The mass of the light particle, /i is much smaller compared to masses of the other mesons 
that all contain heavy quarks, p « m,M,M', so in the following we take fi = 0. For fi ^ 0 the amplitude has a 
slightly different analytical structure. The difference, however, does not affect the behavior near the 2m threshold. 

The contribution of the 2m channel to the discontinuity of the M + fi M' + p 5-wave amplitude, is given by 


A A(s) = B*(s)p(s)C(s)6(s — 4m 2 ). 


( 3 ) 


Here B(s) and C(s) are the 5-wave projections of the amplitudes B(s,t(s, z s )) and C(s,t(s, z s )) describing the 
reactions, M r (jpz)+p(p±) m(qi)+m(q 2 ) and M(pi) +/i(— P 2 ) ^ 7 ,(^ 1 )+^(^ 2 )? respectively. The phase space factor 
p(s) is given by p(s) = yj 1 — 4 m 2 /s. It should be noted that amplitude discontinuity, A A(s) = (A(s+ie)—A(s—ie))/2i 
is in general a complex function of 8 and not equal to Irm4(8), as sometimes assumed [HI ITT) . This has important 
consequences when considering phase of the amplitude, which is often used to discriminate between resonant and 
non-resonant behaviors. 

We reconstruct A(s) from its discontinuity approximated by Eq. © since, in the region of interest, the difference 
between A A(s) and the true discontinuity is expected to be a smooth function of 8, 



ds 


s' — 8 


( 4 ) 


where s tr = 4m 2 . If A(s) is to have a strong dependence on 8 , the numerator under the integral has to vary rapidly 
with 8 . Rapid variations are determined by nearby singularities. As discussed above, if direct channel bound states 
are excluded, the only physically allowed singularities of partial waves are branch points. Thus in order for the integral 
in Eq. Q to develop a singularity the branch points of the numerator have to either appear at the endpoint of the 
integration region or pinch the integration contour. 

We note here, that in PH the amplitudes B(s) and C(s) were chosen proportional to an exponential, exp (—(8 — 
str)/so) with 80 = 0(lGeV 2 ). Such behavior was motivated by a quark model. Non-relativistic quark model 
calculations often involve gaussian wave functions and one could imagine that a calculation, in which the scattering 
amplitude is derived from a diagram involving quark exchanges between interacting mesons, would produce energy 
dependence that falls off exponentially. Such a model, at best, is valid for non-relativistic relative momenta and 
has incorrect crossing-properties that prevents it from being used in a dispersive analysis. A dispersive analysis 
requires that the amplitudes vanish at infinity in all directions in the complex 8 -plane or, in the worst case, grow 
polynomially. Instead, the amplitudes in m have an essential singularity at infinity. Such a singularity is unphysical, 
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it violates causality, which requires the amplitudes to be polynomially bound. If a proper Lorentz covariant quark 
model was constructed the quark exchange mechanism would correspond to ^-channel exchanges that do not overlap 
with exchanges of normal, i.e. quark-antiquark, mesons. Such ” forces” lead to 5-channel partial wave amplitude 
singularities located far to the left from 5-channel unitary cuts and as such do not produce enhancements in the 
5-channel physical region. For example the form exp (—(5 — 5 tr )/5o) might be replaced by (s tr + so)/(s + 5o). The 
latter vanishes in all directions at infinity and with 5o > 0 it is often used to approximate singularities located far 
away to the left from the physical region. When such a form factor is used in Eq. @ (in place of the B and C) it 
results in a normal, not enhanced threshold cusp. 

In the following we show that A(s) develops a pinch singularity through C(s) but not for through B(s). Since C(s) 
is the 5 channel, S-wave projection of the M + fi m + m amplitude 


C(s) = J ^dz s C{s,t(s,z s )) 

we use the t-channel dispersion relation to write, 

1 r°° 

C{s) = ^ d\' 2 Ci(s, A i2 )Q(s, A' 2 ) 

where 


( 5 ) 

(6) 

( 7 ) 


Here C\ is the t-channel discontinuity of C(s,t) at fixed-5 and A^ in is the location of the lowest mass singularity. 
The lowest mass that can be exchanged in the ^-channel is the pseudoscalar or vector (here, for simplicity, assumed 
degenerate and stable) open flavor, e.g. D or D*, meson. In this case 61 ( 5 , A 2 ) = GS( A 2 — m 2 ) (ignoring spin of the 
exchanged particle) and the constant G is given by the product of couplings of the exchanged meson to the (Mm) 
and (m/i) meson pairs. For larger A 2 unitarity in the ^-channel, relates the discontinuity Ci(s, A 2 ) to the product of 
amplitudes M + m m + (n/i) and m + fi —)> m + (nfi) representing the interaction between the external meson 
pairs (Mm) and (m/i), and intermediate states containing, besides the heavy meson, other hadrons, e.g. n light 
mesons, (nil). It will be shown below that the relevant singularities of C(s) originate from a small interval in A 2 that 
corresponds to ^-channel intermediate states containing a few particles. Complicated t-channel intermediate states 
are therefore not relevant when looking for sources of threshold cusp enhancement. 

The key is the function Q(s) defined in Eq. 0. It is the 5-wave projection of the t-channel amplitude describing 
an exchange of an object of mass A. As a function of 5 and the momentum transfer t is given by 


£( 5 , z s ) - M 2 + rn 2 


5 + M 2 

2 


0 


M 2 )V 5 — 4m 2 

2 ^ 




( 8 ) 


Typically, the 5-channel partial wave projection of a £-channel exchange leads to an amplitude with branch points 
located to the left from the 5-channel thresholds and therefore outside the 5-channel physical region. Singularities in 
Q(s ) near the 5-channel physical region can appear, however, because of unequal mass kinematics. Formally Q(s) is 
given in terms of the l = 0 Legendre function of the second type, 


Q(s, A 2 ) 


_ hi* 1 A 2 -t-00 

(M 2 — s)\/4m 2 — 5 A 2 — £+( 5 ) 


(9) 


where t±(s) = £(s, ±1). Determination of Q(s) on the physical sheet where the dispersion relation, Eq. is 
evaluated, corresponds to an appropriate choice of branch cuts of the logarithm and the square roots. 

For A 2 real, the function Q(s) on the physical sheet is obtained by continuing in M 2 from above the real axis, 
M 2 —» M 2 + ie [24]. The function has four branch points, 5 = 0, —00 and 5±(A 2 ), given by the two independent 
solutions of the equation t±(s±( A 2 )) = A 2 . The motion of the s± branch points as a function of A 2 is shown in Fig. [l] 
and Fig. [2] There are four characteristic regions, which we enumerate here, 


i) A 2 > A 2 = (M + m) 2 , 

ii) A 2 > A 2 > A l = (M — m) 2 , 
Hi) A 2 > A 2 > A 2 = M 2 /2 — m 2 , 
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FIG. 1: Motion of the s±(A 2 ) branch points in the complex s-plane as a function of A 2 . The two arrows indicate the limits 
as A 2 —y oo. As A 2 decreases s_ moves to the right above the real axis. For oo > A 2 > A 2 , s+ moves to the right below the 
real axis. For A 2 > A 2 > A 2 , s+ moves to the left below the real axis. As A 2 decreases below A 2 , s+ circles the real axis below 
threshold and moves to the right above the real axis. The point a represents is located near A 2 . Point b is near and above A 2 . 
When M 2 approaches the real axis, s±(A 2 ) pinch the real axis at s = Sb = M 2 m/(M — m). The point c is near and below A 2 . 
For A 2 > A 2 branch points s±(A 2 ) are located on the same side of the real axis and do not produce a singularity in A(s) when 
s approaches real axis from above. Pinching occurs for A 2 in region Hi), i.e. A 2 > A 2 > A 2 . 



FIG. 2: Behavior of the real parts of s±(A 2 ) as a function of A 2 in the region corresponding to the box shown in Fig. [l] 
iv) A 2 > A 2 . 

For large values of A 2 , i.e in region i), Eq. |9| represent a genuine ”force” effect of a particle exchange. Branch 
points of Q(s) are located on the negative 5-axis with cuts extending between the point s = s+(A 2 ) and 5 = 0 and the 
points s — —oo and s = 5_(A 2 ). As illustrated by the arrows in Fig. [lj when A 2 -T oo, 5 + (oo) = 0 and s_(oo) = —oo 
i.e., cuts disappear and the exchange reduces to a point-like interaction. At A 2 = A 2 , the two s± branch points 
collide and as A 2 decreases into region ii) the branch points move into the complex plane. In this region the left 
hand cut extends over the entire negative real axis and along an arc in the complex plane between s+ and s_. As 
A 2 approaches A 2 from above, s+/_ approach the real axis from below and above, respectively. This happens in the 
physical region of the 5-channel. As long as M 2 is slightly above the real axis the two branch points s± never cross 
the real 5-axis. In the limit e —» 0 they pinch it instead. This is the origin of a singularity of A(s). As A 2 decreases 
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c,(s',?i 2 )Q(s;a, 2 ) 



t 


t 



FIG. 3: Left panel: representation of A(s) given by Eq. (10). Right panel: for real s, A(s) is equivalent (up to 

arbitrary normalization), to a Feynman triangle diagram in a scalar theory [27], A(s) = f* da\da2das8(l — a% — ol<i — 
03 ) [ai \ 2 + (ct2 + OL2>)m 2 — ona^/i 2 — 0203 s — 01 O 2 M 2 — ie\ 1 . The internal lines represent scalar particles with masses m 
and A and squares of momenta of the three external lines are M 2 , mu 2 and s, respectively. 


towards A 2 , s + moves to the left, circles the threshold at A 2 = A 2 and, avoiding the unitary cut, when A 2 decreases 
below A 2 into region iv), moves above the real axis. In region iv) the branch point s_ stays above the real axis and 
moves to the right as A 2 decreases. For 5 in the physical region, i.e. approaching the real axis from above we can 
draw the following conclusions. Since the s± branch points never cross the unitary cut, the integral in Eq. © is well 
defined. When A 2 is in region i) or ii) the numerator in Eq. Q has singularities far away from the physical region 
and the only singularity of A(s) is the threshold branch point. For A 2 in region Hi) s and s + are on the opposite sides 
of the integration contour and A(s) develops a pinch singularity of a log-type at s = As A 2 decreases below A 2 
the threshold point, s = 4m 2 is once more the only singularity of A(s). The real axis singularity for A 2 in region Hi) 
can also be traced to the complex singularities of A(s) for other values of A 2 . For example, in region ii) A(s) has a 
singularity on the second Riemann sheet reached by going down the unitary cut from above. As s + approaches the 
real axis from below i.e. A 2 moves from region ii) to Hi), this second sheet singularity moves towards the real axis. 

When M < 2m, which is the case of amplitude B(s), there is no pinching. This is because the condition (M —m) 2 > 
A 2 > M 2 /2 — m 2 requires that the mass of the ^-channel exchange is smaller than the mass of the lightest open flavor 
meson, A < m. Thus B(s) is expected to be a smooth function and in the following we approximate it by a constant. 
In this approximation, A(s) given by 


pOO pC 

A(s) =A 0 ds' / 

J 4ra 2 J X r 


d\ 


/2 


Ci(*',A' 2 V(s')Q(s , ,A /2 ) 


( 10 ) 


This is equivalent to a triangle diagram shown in Fig. [ 3 ] thus Eq. flip] ) has the same singularity structure as the triangle 
diagram in perturbation theory with the branch points s± corresponding to the leading Landau singularities [ 251427 ] . 

To illustrate the points discussed above we consider a model with an exchange of a single particle with mass A 2 . In 
this case C\{s, A /2 ) = 5( A /2 — A 2 ) and, in arbitrary units, we choose M = 10 and m = 4. In Fig. [I] we show modulus 
of A(s) as a function of s for different values of A 2 . The constant Aq in Eq. (10) is chosen such that \A\ is normalized 
to unity at threshold, |A(4m 2 )| = 1. For the choice of masses given above, A^f = 34 and A 2 = 36. For A 2 = 100, i.e. 
in region ii) the numerator in Eq. (10) has branch points in the s-plane far to the left from the threshold and the 
only s-plane singularity of A(s) is the threshold branch point. As A 2 decreases and reaches a point slightly above A 2 
the branch points s±(A) approach the real axis. This produces a bump in s close to s 5 = M 2 m/(M — m) = 8.165 2 , 
which originates from the second sheet singularity discussed above. At A 2 = A 2 pinching begins i.e. in addition to 
the threshold singularity the amplitude develops a singularity at s p = 85 . As A 2 decreases towards A 2 the s + (A) 
branch point moves to the left bringing the pinch singularity closer to threshold. When A 2 reaches A 2 , the singularity 
at s p collides with the threshold singularity producing a sharp peak at s = 4m 2 . When A 2 decreases further the 
pinch singularity moves away from the real axis to the second Riemann sheet reached by moving through the unitary 
cut from below. This singularity is far away from the physical region leaving the point s = 4m 2 as the only visible 
singularity of A(s) in the physical region [28]. Enhancement of the normal threshold cusp and/or appearance of the 
pinch singularity near threshold can happen for A 2 in a limited range. This is important for phenomenology, since it 
implies that only a small part of the t-channel spectral function is relevant. Physical interpretation of the kinematics 
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FIG. 4: Example of |A(s)| obtained using Eq. (10). See the text for interpretation of the results. Color online. 


corresponding to region in) was given in [29] . The strongest effect occurs in the case M = 2m when A 2 = A 2 In this 


case pinch and end-point (threshold) singularities overlap producing effectively a pole of the numerator in Eq. (10) 
at s = 4m 2 . The numerator behaves like a pole only for s > 4m 2 . To the left of s = 4m 2 there is a branch cut of 
infinitesimal length joining the branch points at s±(A 2 ). The physical origin of this pseudo-pole corresponds to the 
limit of the exchanged m particle going on-shell when the fi particle becomes soft. 


III. PHENOMENOLOGY OF THE Z c (3900) AND Z b (10610) 

In the following we apply the formalism described above to the case of the two isovector, Z states, the Z c (3900) 
and the Z&( 10610). The Z c (3900) is observed as a peak in the J/'ipTr mass distribution near the DD* threshold. The 
effect is seen in the decay of Y(4260) to J/'ijjnn. The Z&( 10610) corresponds to a peak in T(nS)n, (n = 1,2,3) near 
the BB* threshold in the reaction T(5S f ) —>• T(lS)7nr. What follows is only an illustration and not a comprehensive 
study. For realistic comparison with the experimental data it would be necessary to include background contribution 
as has been done, for example in M- 


1. Z c (3900) 


In the notation of Sec. [TT] we choose (in units of GeV) M = 4.260, m = (mp + m#*)/2 = 1.936. This gives 
A 2 = 5.33 and A 2 = 5.40, respectively. The t channel exchange of a single D or D* meson corresponds to A 2 = 3.49 
and 4.03 respectively i.e. the single pole in the t-channel gives a function Q(s) with branch points in the region iv) 
and therefore is not expected to significantly enhance the amplitude near the Z c peak. 

At larger values of A 2 the next t-channel singularity of C corresponds to exchange of a Dir or a D*tt system. In 
addition to the continuum these two-body systems have several resonances. The lowest resonance is the Dq( 2400). 
Ignoring the Dir continuum, i.e. using the narrow width approximation, one obtains, 

C 1 Pole = <y(mi,.-A 2 ). (11) 

The Dq exchange falls entirely into the region in), (A 2 = 5.37) and thus produces a sharp peak at s + (5.37) = (3.87) 2 . 
This is shown in Fig. [5] In a realistic case, however, the Dtt channel should by represented by a distribution that 
takes into account the Dn continuum and the resonances. For example, one can approximate the C\ spectral function 
by 

Ci(s, A 2 ) = ImII(A 2 ) (12) 


where 


n (A 2 ) 


m 2 D , 


A 2 - - 

TT 


-f 

Wt t , 


00 dt , p{t')N{t') 


t' — A 2 — ie 


-l 


(13) 

















7 



FIG. 5: Magnitude of the amplitude A(s) for T(4260) —> JNtktv. The dashed line corresponds to the narrow width 

approximation of Eq. (11) and the solid line corresponds to Eq. (12). The amplitude and the data are normalized to unity at 
threshold 



FIG. 6: Comparison between C\ given by Eq. (12) with the imaginary part of the Breit-Wigner formula, oc [m 2 — t — mV] 
with PDG values for mass m — 2.318 and width T — 0.267 of the Dq. Normalization of the amplitudes is arbitrary. 


J12 


k reproduces 
width of the 


tt r = {tyid 0 + 777/7T) 2 . Choosing friD* = 2.67 and N(t) = gA. 2 /(t + A 2 ) with g = 1.5, A = 10, Eq. 
the Breit-Wigner line shape with parameters corresponding to the physical values for the mass and t 
Dq [30] . The comparison between Ft (t) and the Breit-Wigner formula is shown in Fig. [6] The reason why C\ has no s- 
dependence in this model is because we ignored spin of the t-channel Dir system. This is justifiable since s-dependence 
from the spin of an exchanged object is a smooth function. Integration over the ^-channel mass distribution of the Dtt 
spectrum, as expected, has the effect of smearing the sharp peak obtained in the narrow resonance approximation. 
Nevertheless the cusp at threshold remains enhanced by presence of the nearby branch points Q(s). There are other 
D-mesons that can contribute in the mass region of the Z c (3900) and in a phenomenological analysis of the data 
should be considered. For example the Di(2420) is much more narrow that the Do (2400), however, its nominal mass 
places it in region ii). 

In Fig. [7] we compare the real and the imaginary part of A(s) obtained in the narrow width approximation, (c/. 
Eq. ([IT])) with those computed using Eq. (12). The amplitudes are different from the case when the left hand cuts 
are far away from the threshold branch points (c/. Fig. |8|. The numerator under the integral in Eq. @ is a complex 
function of s' and so is the amplitude below threshold, s < 4m 2 . Above threshold, in the narrow width approximation, 
for A 2 in region in), the branch points of the numerator are located on the opposite side of the real axis producing an 
imaginary amplitude for s > A 2 . 













2 



| - Re(Pole(t)) 

1 

1 

j -lm(Pole(t)) 

1 

1 

\ - RetC^t)) 

1 

1 

1 

1 

1 

1 

j - ImfC^t)) 

1 

h 

_—"— J/% 

\ 

- /' 

/ i 
/ 

_—— / 


-^ 

_ 


0 -- 1 --‘- 

3.8 3.9 4 

s 1/2 


FIG. 7: Real and imaginary part of A(s) computed using Eqs. (labeled pole) and using Eq. |l2| ), labeled C\. 



FIG. 8: Modulus (solid black), real (solid red) and imaginary (dashed red) parts of the amplitude for Y (5 S) —> T(1S)tttt 
obtained using the formula of Eq. (12) with fh D * fh B * » 6.036 and g — 7.5. As discussed in the text, the relevant range of 
A 2 corresponds to the region ii) which produces a singularity in the amplitude located to the left of the unitary branch point. 
The real and imaginary parts thus exhibit behavior that is typical for an amplitude with a left hand side singularity. 


2 . Z b (10610) 


We apply the same analysis to the decay T(5S) -T T(1s)7T7t. In this case we use M = 10.86, and m = (tub + 
rriB*)/2 = 5.302. The t channel exchange of B or B* corresponds to A 2 = 27.87 and A 2 = 28.36, respectively, which 
are below A 2 = 30.86. That is, just as in the case of the Z c (3900), one-particle-exchange contributions fall into the 
region iv) where threshold cusp is not enhanced. In the relevant mass range, the Bit or B*tt spectrum have only 
one known resonance, the B^ with mass m = 5.698 and width T = 0.128. In the narrow width approximation this 
corresponds to A 2 = 32.47 and is above the pinch region, which corresponds to a narrow mass window between, 
30.86 < A 2 < A 2 = 30.89. It is only the low mass tail of the Bit distribution that overlaps with this very small region 
where the pinch singularity of A(s) occurs. For example, the mass of the Bit system where pinching overlaps with 
threshold corresponds to A = A c = 5.555, which is less than one resonance width, T, below the Breit-Wigner mass of 
the B*j resonance. Therefore, even though the resonance is away from the pinch region it can nevertheless enhance 
the normal threshold cusp. This is illustrated in Fig. [8j The comparison between Eq. (12) with Bit system parameters 
and the corresponding Breit-Wigner formula is shown in Fig. [9] 
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FIG. 9: Comparison of the formula in Eq. 


obtained using the PDG parameters m B * = 5.698 


(12) for the bottomonium with the imaginary part of the Breit-Wigner formula 
= 0.128. 


r B ; 


A. Summary 

Motivated by the recent discovery of the XYZ peaks in charmonium and bottomonium spectra and by theoretical 
models that propose to explain their origin in terms of normal threshold cusp enhancements, we preformed a systematic 
analysis of 2-to-2 amplitudes in the kinematic region of interest. We found that amplitudes containing quarkonium 
with mass above the open flavor threshold, the amplitude C in notation of Sec. [TT| can have singularities in the 5-plane 
that enhance the threshold cusp. For the enhancement to occur, the absorptive part of the amplitude has to have 
a singularity on the second Reimann sheet close to the physical region. When this happens the dispersive integral 
develops a pinch singularity and when the singularity lying below the real axis coincides with threshold, it enhances 
the normal threshold cusp. This, however, only occurs provided the ^-channel spectral function is large in a narrow 
window of masses corresponding to the region Hi). An amplitude with the quarkonium mass below the open flavor 
meson (i.e. the amplitude B ) is not expected to be significantly enhanced. This can be understood by considering, for 
example, the case M = 2m as discussed in Sec. [TIJ We thus find that in the kinematical region studied here, threshold 
cusps are enhanced by the same type of left hand cut singularities as present in triangle graphs. Such a mechanism 
has been proposed in [11]. In Ref. [17] both the B and the C amplitudes were assumed to be enhanced, contrary 
to what is expected based on the arguments presented here. Furthermore, in the the previous studies the amplitudes 
were given an analytical form that is not physical. In relation to the XYZ physics, triangle diagrams were used in 
H2HH. These works did not, however, present a systematic study to explain which exchanges are relevant. In fact, 
it was assumed that it is the exchanges in the amplitude B and not C that are important. For example in m an 
exchange corresponding to our amplitude B was found to produce an effect that could in fact describe the data quite 
well. It would be interesting to analyze the singularity structure of that amplitude and compare with our predictions. 
We have performed a preliminary analysis of J/7r and T(1S)tt amplitudes and shown that the corresponding triangle 
singularities can potentially produce enhancement in the amplitude qualitatively consistent with the data. It should 
be noted, however, that such a mechanism cannot be considered as the final answer. If singularities of C were far 
away to the left from the physical region they would represent a genuine ”force” i.e. virtual particle exchange, an 
approximation to the left hand cut that can be unitized in a standard way m- This point was also raised in m 
against findings obtained with the model of PH- The situation is however, even more complex since singularities in 
the physical region found in the triangle diagrams are due to opening of new channels and these will have to be taken 
into account in the unitary relation. 
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